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Abstract 

In the present work we deal with the dynamics of wave packets in a two-dimensional crystal 
under the action of magnetic and electric fields. The magnetic field is perpendicular to the plane 
and the electric field is on the plane. In the simulations we considered a symmetric gauge for the 
vector potential while the initial wave packet was assumed to have a gaussian structure with given 
velocities. The parameters that control the kind of time evolution of the packets are: the width 
of the gaussian, its velocity, and, the intensity and direction of the electric field as well as the 
magnitude of the magnetic field. In order to characterize the kind of propagation we evaluated 
the mean-square displacement (MSD), the participation function and which is more important 
we were able to follow the wave at different times, which allowed us to see the time evolution 
of the centroid of the wave packets. A novel effect was observed, namely, the dynamics is such 
that the wave function splits into two or more components and reconstructs successively as time 
goes, vortices are forming. To our understanding this is the first time such an effect is reported. 
As for the inclusion of the electric field, we observe a complex behavior of the wave packet as 
well as note that the vortices propagate in direction perpendicular to the applied electric field. A 
similar behavior presented by the classical treatment, In our case we give a quantum mechanics 
explanation for that. 

PACS numbers: 73.20.Jc; 72.20.My 73.23.-b 



1 



I. INTRODUCTION 



In the present work we deal with the problem of the behavior of wave packets in a two 
dimensional square lattice under the action of magnetic and electric fields. The magnetic 
field is perpendicular to the lattice, while the electric field is in the plane. 

We can mention the (pioneering) works done on the subject of wave propagation in 
low dimensional systems that have attracted the interest, since the early days of quantum 
mechanics.-^ 1 ^ 

The subject of carriers in a two dimensional structure under the action of external mag- 
netic and electric fields, has aroused intense interest since it has become recently experimen- 
tally accessible .-^^ 

We have found very interesting properties of the time evolution of initial wave packets 
that were assumed of a gaussian structure with a given velocity. We analyzed gaussians 
with different dispersions which in turn determine the type of propagation a wave packet 
will present. Another parameter which has a direct influence on the wave packet behavior 
is the assumed initial velocity. Obviously, the magnitude as well as the direction of the 
electric field has also a direct influence on propagation. Clearly, one notices the wavy nature 
of the solution of the time dependent Schrodinger equation but, a novel effect is observed, 
namely, the successive splitting and reconstruction of the wave function in two or more 
components as time goes. We observe the formation of vortices due to the joint effect of the 
crystal potential and the external fields. To our understanding this novel effect has not been 
reported until now. In order to comprehend the characteristic of propagation we resort to 
the study of the trajectories in reciprocal space since they are connected with the ones in 
direct (coordinate) space by a rotation of 7r/2. 

A pioneering work dealing with the motion of an electron in a 2D lattice potential su- 
perimposed to a magnetic field is due to Peierls^, that considered an effective single band 
Hamiltonian arising from a tight-binding dispersion relation. As a result of this model the 
single Bloch band is split into magnetic subbands according to the number of flux quanta 
that pierces the unit cell of the 2D lattice. At the same time, the parameter a = $/$ being 
the ratio between the magnetic flux through the unit cell to the quantum of flux, controls 
the kind of propagation of wave packets in a lattice under a uniform magnetic field. For a 
rational value we recover translational symmetry, with an enlarged unit cell which in turn 
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makes possible for a packet to propagate in the sample. On the contrary for an irrational 
value of a we face the problem of incommensurability of the potential which in turn produces 
a localization of a wave packet in a definite region of the lattice^. Recently, the electronic 
spectrum of a two-dimensional quantum dots array under magnetic and electric fields was 
presented^. 



II. THE MODEL 

The action of a magnetic field is analyzed along the Peierls model^, which consists in 
taking a dispersion relation for a square lattice 

E(k) = 2W (cos k x a + cos k y a) (1) 

and replacing in it the quasimomentum k by 

hk =^ -ihV - eA/c (2) 

to obtain a model Hamiltonian. In the present work we used the symmetric gauge for the 
vector potential: 

A=|(-2/i + xj + 0k) (3) 

The classical paper of Hofstadter 12 showed that such an approach leads to a spectrum 
as function of the magnetic field that presents a fractal structure, the so called Hofstadter 
butterfly. On the other hand, Hall measurements on GaAs/AlGaAs superlattice provided 
evidence of the existence of the structure of Hofstadter's butterfly^. Since this model is of 
a single band, is limited to analyze systems of large gaps and/or magnetic intensities such 
that no interband transitions occur. Since we want to study the kind of propagation of a 
particle in such a system, we expand the wave function in the Wannier representation: 

mn 

where the ket | mn) is the ket associated with the corresponding site^. Next, we assume a 
discrete set of coordinates such that x = ma, and y = na. The time dependent Schrodinger 
equation in the Wannier representation becomes the set of equations: 

dg mn _ , inan , -inan , -inam , inam \ 

dt ~ K9m+l,n e i 9m-l,n e + 9m,n+l 6 -|~ y m ,n-l e ) 

+ 9mn (e mn + eE x am + e E y an) (5) 
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where a = $/$o is the ratio between the flux through the unit cell in the (x, y) plane to the 
quantum of flux 10 $ = he/ e, e mn are the on-site energies and W is the hopping term. 

We have used the Runge Kutta method of fourth order to integrate the equations of 
motion. In order to solve the time dependent Schrodinger equation we chose as initial 
condition a gaussian wave packet with a certain width and a given velocity: 



( x,y | *ff(t = 0) ) = expi (k • r) — = exp 

awn 



(x - x ) 2 + {y- y f 



(6) 
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We decided for this since it is more realistic to assume that the injected electron is not 
extremely localized. One difficulty faced during the calculations was to decide the right size 
of the lattice in order to avoid boundary effects. To be specific, we have taken: a lattice of 
400 x 400 sites, a = 0.0242 which corresponds to a magnetic field of intensity B = IT, the 
dimensionless units of time r = Wt/h, corresponds to 0.026ps, the dimensionless units of 
electric field are W/ea, such a unit is equivalent to 12.5 kV/cm, and all lengths are in units 
of the lattice parameter a. After solving the set of equations we constructed the following: 

i) the mean square displacement (MSD) 

<r 2 >(t) = \9mn(t) | 2 (nV + mV) (7) 

mn 

ii) we follow the centroid of the wave packet by evaluating the following quantity, which 
give us the amount of the displacement form the initial position of the particle: 

< Ax > (t) = J2 ( U ~ n o) Yl I 9mn{t) | 2 (8) 
n m 

<Ay>(t)= ( m - m o) Yl I 9mn(t) | 2 (9) 

m n 

iii) the participation function^ 



P(t) 



^ ^ I 9mn(t) 



(10) 



An interesting feature of this function is that it indicates the sites that participate in the 
wave packet. At the same time it presents an abrupt decline once the packet reaches the 
boundary of the lattice, in this way we can note the presence of size effects. We followed 
the Anderson^ criterium for analyzing diffusion, namely, we can conclude that diffusion has 
occurred if at t — > oo the Wannier amplitudes at the starting sites go to zero. If, these 



amplitudes remain finite decreasing rapidly with distance, we say we have a localized state. 
We also plot the wave packet as it evolves in time which tells us the kind of propagation for 
the different cases in study. More than that, by looking at the displacements of the maxima 
of the packet we can infer the kind of trajectory a particle would describe. Besides that, 
we follow the time evolution of the centroid of the wave packet which gives complementary 
information. 

III. THE SPLITTING OF THE WAVE PACKET: VORTICES FORMATION 

We would like at this point to describe a surprising behavior observed during the time 
evolution of wave packets. First, we note that their gaussian structure means that in recip- 
rocal space we have a certain dispersion in k, which implies that several wave vectors will 
participate in the evolution of the wave function. This plays an important role as long as we 
are considering the wave vector associated with the velocity of the initial wave packet, lying 
along and around the lines of zero energy in the Brillouin zone of the square lattice. In Fig. 
1 we show the lines of constant energy where the arrows signal the orbit in reciprocal space 
described by the wave vector. 

As it is well described in the books of Solid Stated, the quasimomentum satisfies an 
equation of motion analogous to the classical one; where one takes the group velocity v = 
in the expression of the Lorentz force, which in turn determines that the wave vector moves 
along the lines of constant energy. First, consider k inside the region of energy zero, but 
close to its boundaries, it will describe a clockwise orbit in reciprocal space. As for k outside 
but close to the line, we get another trajectory described in counterclockwise sense. This 
will result in the appearance of vortices rotating in opposite directions when describing the 
evolution of the wave packet in direct space, where the "orbits" are obtained after rotation 
by tt/2. 

What is very interesting is to consider k on one of the lines of zero energy. In such a case 
due to the dispersion in k because of the gaussian structure of the initial packet, we discussed 
above, we have to take into account k values inside and outside the line, symmetrically 
distributed. Consequently, part of the wave (half of it) will describe a trajectory in the 
clockwise sense and the other half of the wave in the opposite sense, the resulting movements 
is similar to a swimmer doing breaststroke. As time goes, the wave packet starts moving 
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in the direction of the velocity and then splits into two components that will join in a 
single component and so on. Consequently, two vortices with opposite angular velocities are 
formed. The wave remains stationary in a certain region of the lattice. See Fig. 2 where 
this remarkable effect is shown. Assuming now k with components (ir/a,0), as shown in 
Fig. 1, four squares in reciprocal space participate in the movement of the packet reinforcing 
each other such that at a certain time the wave appears like is shown in Fig. 3. To get the 
evolution of the wave packet one has to follow the arrows in each of the four regions. Again, 
we observe the splitting and reconstruction of the wave, this time in several components. In 
this case, several vortices are present. 

One encounters a similar effect when considering a 1-D crystal under the action of a dc 
electric field. When considering as an initial state a well localized state at a site in the 
system, the evolution of the packet is such that it is split into two symmetric parts, which 
oscillate with the Bloch frequency in opposite directions^. 

Let us assume k inside the region limited by the lines of zero energy but close to one 
of them. The wave function will split into two components in such a way that one part 
of the wave, the major part, will rotate clockwise while the rest will do counterclockwise. 
As for k outside the region but close to one of the lines, the reciprocal is true, the major 
part of the wave will perform a rotation counterclockwise. To illustrate this effect, consider 
now the wave vector of the initial wave packet, inside the region of the lines for e = but 
close to one of it, we took for example ka = (1.4, 1.4). As said above, the wave is split in 
two components, where the bigger part rotates clockwise and the smaller part rotates in the 
opposite sense. In this case this case it results in the appearance of two asymmetric vortices. 
See Fig. 4. 

IV. THE EFFECT OF A DC ELECTRIC FIELD 

We consider now the inclusion of a dc electric field in the equation of motion for the 
Wannier amplitudes. As a general trend we observe that the packet will propagate in 
a more complex way, but always in direction perpendicular to the electric field, as it was 
shown previously^. This behavior is also present in the classical treatment of the problem^. 
From the view point of Quantum Mechanics, we understand this behavior since the electric 
field breaks the degeneracy of the on-site energies along the direction of the applied field, 
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inhibiting hopping between these sites. First we take k on one of the lines of zero energy, 
for example (ir/2a, n/2a) and the electric field along the diagonal, E = (0.1,0.1), while the 
initial packet has a — 1. The wave is split while propagating but, due to the presence of 
the electric field, one part of the wave proceeds with a greater velocity, with a centroid 
trajectory similar to the one a snake performs when moving around. The other part that 
moves in opposite direction, remains close to the starting point. See Fig. 5 where we also 
show the centroid trajectory, and the MSD and participation as functions of time. 

For the same configuration but for a = 2, we observe that the wave is not split and follows 
a trajectory quite similar to the classical one, i.e., the more so, the more extended is the 
initial wave packet, i.e., the greater is sigma. This comes about since a greater dispersion in 
direct space is related to a smaller one in reciprocal space. See Fig. 6. For the case a = 3 we 
confirm this behavior of the propagation of the wave. It is interesting to mention that the 
trajectories for a = 2 and 3 are exactly the same. By increasing the intensity of the electric 
field we obtain a displacement of the centroid along a trajectory similar to the former one, 
the only difference is that in the case of stronger field the amplitude of the displacement as 
well as the period of the oscillations are reduced. This field effect is shown at the bottom 
left of Fig. 6. 

A very peculiar effect is obtained by considering the electric field with components 
(—0.1,0.1) and taken k = (n/2a, n/2a). In this case we note that the "swimmer" by doing 
breaststroke can displace itself along the perpendicular direction to the applied electric field. 
The vortices displacement is such that the centroid moves along a strait line, as shown in 
Fig. 7. This should be compared with the case without electric field in which the packet 
remains stationary as shown in Fig. 2. 

V. CONCLUSIONS 

We show in this work a novel effect of wave packet propagation in a two dimensional 
crystalline system under the action of combined magnetic and electric fields. Using the 
Runge Kutta method of forth order, we integrate the equations of motion in the Wannier 
representation, assuming as initial condition a gaussian wave packet with given velocity. 
We found the very interesting behavior of the wave function, namely, that by taking the 
initial velocity with the associated wave vector close to the lines of zero energy, the wave is 
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split and reconstructed as time goes, with the appearance of a series of vortices. This effect 
comes about since we used a gaussian as an initial condition. This in turn, implies that 
one has to take into account a dispersion in the reciprocal space, so there are contributions 
of the lines of constant energy on both sides of the zero energy line, as explained above. 
Without the presence of the electric field, the wave remains in a definite region of the 
lattice. The inclusion of the electric field, produces a displacement of the vortices along 
the perpendicular direction of the applied field, showing a more complex behavior since 
now, besides the displacement, the wave is split and reconstructed as time goes. As for the 
centroid trajectory, for a configuration such that the applied electron field is parallel to the 
initial velocity, it describes a trajectory similar to the way a snake performs when moving 
around. In the classical treatment, the trajectory is a trochoid. For other configurations, i.e., 
for the wave packet with initial velocity perpendicular to the field, the centroid trajectory 
is a strait line. One last comment deserves to be made, the effects we described are the 
results of taking into account that the particle, besides being under the action of the fields, 
is subjected to a 2D crystal potential as well. 
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FIG. 1: We show the lines of constant energy for the square lattice. The solid lines are the ones 
corresponding to zero energy. The arrows indicate the orbits in reciprocal space. 
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FIG. 2: The time evolution of the wave packet for the following parameters: a = 0.0242, dispersion 
a = 2, and ka = (vr/2,7r/2). 
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FIG. 3: The same as Fig. 2 but for ka = (tt,0). 
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FIG. 4: The same as Fig. 3 but for ka = (1.4, 1.4). 
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FIG. 5: The time evolution of the wave packet for the following parameters: E = (0.1, 0.1), a = 1 
and ka = (ir/2,ir/2). On the left, MSD and participation as functions of time, and the centroid 
trajectory which resembles a trochoid. 
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FIG. 6: The same as Fig. 5 but for a = 2. The bottom left shows the trajectories of the centroids 
for two field intensities: E = (0.1,0.1), and E = (0.2,0.2). Note the different periods when varying 
the field intensity. See the text. 
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FIG. 7: The same as Fig. 6 but for E = (—0.1, 0.1). Note the strait line described by the centroid. 
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